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Relativistic  Broadening  Near  Cyclotron  Resonance 

Kaya  Imre    and  Harold  Weitzner 

Courant  Institute  of  Mathematical  Sciences 
New  York  University,  New  York,  New  York  10012 

ABSTRACT 
Relativistic  broadening  of  absorption  (or  emission)  lines  near 
cyclotron  resonance  in  a  warm  plasma  is  investigated  using  the 
linearized  relativistic  Vlasov-Maxwell  system.   The  unperturbed 
state  is  assumed  to  be  isotropic,  but  not  necessarily  Maxwellian. 
The  expansion  parameter  is  n  =  v  /c,  v   being  the  electron  thermal 
speed.   It  is  assumed  that  the  wave  frequency,  plasma  frequency,  and 
cyclotron  frequency  are  all  comparable  in  magnitude,  and  the  refrac- 
tive index  n  =  0(1).   The  parameter  a  =  n/n   is  of  arbitrary  order, 
thus  the  results  are  uniformly  valid  for  all  values  of  oblique  pro- 
pagation angles,  although  the  relativistic  effects  are  negligible 
for  a<<l.   The  dielectric  tensor  is  reordered,  and  the  dispersion 
relation  appropriate  for  this  problem  is  derived  to  the  lowest  sig- 
nificant order  in  n •   The  results  are  expressed  in  terms  of  the 
readily  calculable  (generalized)  plasma  dispersion  function  Z.   In 
the  Maxwellian  case  the  results  are  algebraic  in  Z,  and  unlike  the 
previously  published  results,  they  do  not  involve  infinite  integrals 
or  series  imposed  on  Z,  thus  leading  to  simple  and  efficient  evalua- 
tions.  The  case  of  perpendicular  propagation  is  obtained  by  taking 
the  large  a    limit.   Some  inconsistencies  in  the  literature  dealing 
with  the  extraordinary  mode  are  resolved. 

PACS  Index:   52.25.Ps,  52.60.+h 
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I.  INTRODUCTION 

The  importance  of  the  relativistic  effects  on  the  broadening  of 
the  cyclotron  emission  lines  has  long  been  recognized.   Trubnikov 
was  first  to  point  out  as  early  as  1956  that  the  Doppler  broadening 
of  spectral  lines  emitted  by  a  magnetized  plasma  in  the  direction 
perpendicular  to  the  field  is  a  purely  relativistic  effect,  which 
must  be  taken  into  account  even  for  relatively  low  temperature 
plasmas.   Although  he  did  not  investigate  this  phenomenon  in  his 
well-known  paper  ,  he  derived  two  distinct  expressions  for  the  plasma 
dielectric  tensor  based  on  the  relativistic  linearized  Vlasov-Maxwell 
system  of  equations.   The  first,  involving  an  infinite  series  in 
terms  of  the  Bessel  functions  J  ,  appears  to  be  particularly  useful 
when  the  spectral  lines  are  well  resolved  (low  temperatures  and  low 
harmonics),  whereas  the  second  form,  which  is  expressed  as  an  infinite 
integral  containing  the  modified  Bessel  functions  K  ,  is  preferable 
when  the  spectrum  is  continuous  (high  temperatures  and  high  harmonics), 
as  was  stated  in  Ref .  1  .   The  second  form  has  the  shortcoming  that 
the  unperturbed  state  has  already  been  assumed  to  be  the  relativistic 
Kaxwellian  distribution  so  that  the  velocity  space  integrations 
could  be  carried  out.   No  such  assumption  need  be  imposed  in  the 
derivation  of  the  former  representation. 

Various  heating  schemes  subsequently  proposed  for  plasmas,  which 
utilize  the  cyclotron  resonance  process,  have  initiated  further 
research  in  which  the  energy  absorption  mechanism  associated  with 

relativistic  broadening  has  been  investigated  extensively  in  certain 

2 
ranges  of  the  plasma  parameters.   Dnestrovskii  et  al.   attacked  the 
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relatively  simpler  problem  of  perpendicular  propagation  in  a  weakly 

3 
relativistic  plasma,  and  Shkarofsky-^  extended  their  results  to  the 

case  of  nearly  perpendicular  propagation  where  the  relativistic 
effects  are  still  significant.   In  both  of  these  studies,  as  well  as 
many  others  which  followed  them  ,  the  integral  representation  of  the 
dielectric  tensor  was  utilized,  thus  restricting  the  system  to  be 
near  thermal  equilibrium.   The  main  purpose  of  the  present  work  is 
to  remove  this  restriction  by  reformulating  the  problem  starting 
from  the  infinite  series  representation  for  the  dielectric  tensor 
for  weakly  relativistic  plasmas.   For  simplicity,  the  unperturbed 
state  is  assumed  to  be  isotropic,  but  not  necessarily  Maxwellian, 
and  only  the  electron  cyclotron  resonance  is  considered.   The 
extension  of  the  present  method  to  the  case  of  ion  resonance,  amd 
to  the  higher  harmonics,  is  straightforward.   As  was  done  in  pre- 
vious studies,  the  plasma  is  assumed  to  be  warm,  i.e.  only  the 
leading  contributions  of  the  parameter  nr^   =  v  /c  are  considered,  where 
V   is  the  electron  thermal  speed  (2T  /m  )^,  and  c  is  the  speed  of 
light.   The  electron  Larmor  radius  is  also  assumed  to  be  small 
compared  with  the  wave  length,  and  their  ratio  is  assumed  to  be  of 
first  order  in  Tl.   The  refractive  index  n  =  kc/co  =  0(1). 

The  formulation  presented  is  uniform  in  terms  of  the  parameter 

0^=  vy(n„c),  (1) 

n„  being  the  longitudinal  refractive  index.   The  case  of  perpendicular 
propagation  (oC  >cd),  and  the  nonrelativistic  limit  (  ^^<  1)  are 
obtained  by  taking  appropriate  limits.   Section  II  contains  these 
results,  as  well  as  the  evaluation  and  ordering  of  the  dielectric 


-Zj.. 


tensor  near  the  fundamental  electron  resonance,  viz.,  1  -\n\/a)  = 

2 
0(T|  ),  where  il  denotes  the  electron  cyclotron  frequency. 

In  Sec. Ill,  the  zero  order  dispersion  relation  is  derived.   This 
relation  differs  from  the  cold  plasma  dispersion  relation  by  the 
presence  of  a  function  A  which  generates  the  relativistic  Doppler 
effects.   Section  IV  is  devoted  to  the  special  case  of  the  Maxwellian 
unperturbed  state.   It  is  shown  that  in  this  case  the  results  can  be 
expressed  algebraically  in  terms  of  the  standard  dispersion  function 
Z,  thus  providing  a  considerable  improvement  over  the  form  given  in 
Ref.  3i  which  contains  an  infinite  complex  integral  or  infinite  series 
involving  the  function  Z.   Equivalency  of  the  present  results  to  those 
previously  given  is  demonstrated  to  order  T|  .  The  study  of  the  per- 
pendicular extraordinary  mode,  which  can  be  found  in  Sec.  V,  requires 
a  second  order  analysis  for  the  investigation  of  the  lowest  order 
Doppler  effects.   Certain  discrepancies,  stemming  from  an  improper 
evaluation  of  the  second  order  terms,  are  found  in  the  previously 
published  results.   Section  V  also  contains  a  study  of  the  nearly 
perpendicular  ordinary  mode  based  on  the  results  presented  in  this 
work.   In  the  part  C  of  the  same  section,  the  Doppler  broadening  of 
the  Nth  electron  cyclotron  resonance  is  studied,  and  particular 
attention  is  given  to  the  second  harmonic  extraordinary  mode. 
A  summary  of  results,  and  a  final  discussion  are  found  in  Sec.  VI. 
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II.  EVALUATION  OF  THE  DIELECTRIC  TENSOR 


The  dielectric  tensor  for  a  plasma,  which  can  be  expressed  for  a 
plane  wave  of  the  form  exp(ik«x  -  iyt)  as 

6=1^  (4Ki/u))  0-,  (2) 


0-  being  the  conductivity  tensor,  can  be  derived  from  the  linearized 

relativistic  Vlasov  equation,  as  was  done  in  many  texts  on  plasmas. 

It  is  convenient  to  adopt  the  (■*",-,(!)  coordinate  reference  system, 

which  is  defined  by  means  of  the  transformation  matrix  P,  and  its 

inverse  P~  ,  where 

0 
P  =   1  -i   0  ,   and 


1  i 
1  -i 
0   0 


,-1  - 


z         z 

-ii  ii 


0 
0 

1 


(3) 


0    0 

representing  contravariant  and  covariant  transformation  rules  applied 
to  the  cartesian  system  in  which  the  z  axis  is  selected  along  the 
external  magnetic  field.   One  has,  in  this  representation, 


where 


T 


PCOD 


(^«-P   • 


.2  /  2 


u  =  v-^,  y=  (i^uVc^)s  a;  =  (cj/j2J|f  -  k  u  /a 


PJ     ""J  j'""J'   J    "j"0""j"'   ^0   "^QX 


OJ^:;  =  4-7Cn:;eym,,  Q^  =  e.B^/m^jC,   B^  =  B^'S  , 


(^) 


m.,  e.,  n.  denoting  the  rest  mass,  electric  charge,  number  density 

J       J       J 

(respectively)    of  the    jth  species,    and  [  with   \d-^uF(u)  =   1    J 


iu  ^  j2 


i„  bF     7.M^  ^        . 
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In  Eq.(5)i  the  unperturbed  distribution  function  is  assumed  to  be 
isotropic,  so  that  Uj_bF/dU|  =  U||bF/du_Lf  and  the  argument  of  the  Bessel 
functions  J  is  kj_Uj_/ll..   The  azimuth  of  the  propagation  vector  k_ 
is  denoted  by  ^  .   It  is  convenient  to  normalize  F  so  that  F(u)d-^u -> 
F(w)d-^w,  where  w  =  u/v . ,  v.  denoting  the  thermal  speed  of  the  jth 
species.   Near  the  electron  fundamental  resonajice,  one  may  neglect 
the  ion  dynamics,  and  introduce  the  parameter 

I   =  (aj-1IlJ)A„v^  =  (1  -|ii^l/cj)(K7]-2,  (6) 


so  that 


^^    0)  ,   cyw^ 


indicating  that  the  resonance  condition  occurs  at  p=-l.   One 
also  assumes  that  the  electron  Larmor  radius  is  small  compared 
with  the  wavelength,  so  that  the  Bessel  functions'  can  be  replaced  by 
the  first  few  terms  of  their  power  series  representations.  With 
^j.'^p^e  ~  ^^''l^'  °^^  obtains 

where  p  ?    ? 

D  =  Wjj  -  ^  -  iofw''.   and   X  =^pe/-^e  * 
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2 
The  7]   terms  are  kept  m  C^.^.i  since  they  are  needed  in  studying  the 

broadening  of  the  perpendicular  extraordinary  mode  in  Sec.  V.   It  is 

clear  from  Eq.  (8)  that  €__  =  OC-q"^) ,  €_,(  =  OCif)"^),  and  £^^,  €  ^_,  ^^^^^ 

are  all  of  zero  order  in  to.   Hence,  by  writing 

e=  71-2^^-2)  ^  -n-M-l)  .  e(0)  ^77^(1)  .  ...,  (9) 

one  obtains  the  nonzero  elements  up  to  first  order  in  ti  as 
6^-2)  =  Xa 


6V^eif=2  6V.^)e-iy'=Xn.a,,,. 


ei°)  =  1  -  ix. 


(0)     = 


(10) 


1   -   Xnia^^Q-  (Zf/oOXaQ^Q, 

ej,°)    =    1    -   X   -   iXniaQ^2' 
^(0),-2i^.,(0)^2i^.    iXnia^.O' 


where 

-n,m   "  "^^Zn.m   ""  2^^  Hf  ^2n+2,m'  '2^)  V^-^"   '  "^H  ^x 


4.     0^2       d      „  a       <<     f  d\    2n+l  ,  m  ^    , 

a„   _    =  (\^.„   _        27rr2   d^  ^2n+2,m-    ■2n+2  J     D  ""-«• 

(11) 


Z^  ^    =   fd^w^w";  F(w)/D. 
n,m       J  -^    u 

It  is  shown  in  Appendix  I  that  Z„^  ^  and  its  derivative  (thus  a   ) 
can  be  expressed  in  terms  of  Zq  ^,    giving 

^0,0  =^^0.0  -  °^\.l' 

^0.1  =  -^-'^^^O.O  *  20tZo^i  -  i30.%^2.  (12) 

^0.2  =-2«^^0.1  "  30^^0.2  -  2'^%.3- 
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Furthermore ,   the  functions  Z_     (thus  Z    )  can  be  expressed  in 

0 ,m         n ,m 

terms  of  the  integrals   (see  Appendix  I  for  details) 

i^  =  [dCi;2^"^z(r;).  (13) 

c 

where 

Zit)    =  7l/ ^  dwF(w)/(w-  ^)  (14) 

with  F(w)  extended  as  an  even  function  for  w<^0,  i.e.  F(w)  =F(-w). 

The  normalization  of  Z{1^)    is  chosen  so  that  when  F  is  I.laxwellian, 

-3/2   -w^ 
namely  F  =  Tt  "^   e    ,  it  reduces  to  the  standard  plasma  dispersion 

function.   This  case  will  be  further  studied  in  Sec. IV.   One  obtains, 

ith  <g;.=  Jd^wgCw)  x<^(w), 


w 


^0,0  =  -  |<^"^>  -  2Io' 
0,1     ^^d  /SO     1  (^^^ 


z 


0,2 


4<w-2>(^-f)-2|%-2o<fl,-K%. 


3o<.  <^ 

The  following  comments  regarding  Eqs.(13)  and  (l4)  are  in  order. 
Although  the  integral  in  (l4)  defines  a  sectionally  holomorphic 
function  in  the  complex  t,   plane  (cut  through  the  real  axis),  it  is 
customary  to  consider  only  the  upper-half  branch  and  analytically 
continue  it  to  the  lower-half  plane,  since  the  lower  branch  can  be 
obtained,  if  needed,  by  Z*(2;*),  which  follows  from  the 
Schwartz  reflection  principle.     Moreover,  since  F(w)  is  an  even 
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f unction  of  w,  the  function  defined  by  this  integral  is  odd  in  t- 
Therefore,  the  values  of  Z(0  ^^   "the  first  quadrant  are  related  to 
those  in  the  second  quadrant  by  Z(^)  =-Z*(-^*). 

The  contour  C  in  Eq.(13)  is  defined  as  the  image  of  the  interval 
\^|  ^  1  under  the  mapping  o/^  =  ^  "^  (jJi  -  (J^)^,  where  jj^  =  2f(y  ,  (see 
Appendix  I).   The  contour  lies  on  the  real  C  axis,  except  when  ]  ^  j  ^^fi^,  i 
which  case  G  lies  on  the  upper  half  of  the  circle  of  radius  {2^/<x)'^i.  R^. 
If  ^ is  allowed  to  approach  the  real  axis  from  below,  then  the  entire 
path  C  remains  in  the  upper  half  plane,  which  enables  one  to  use  the 
conventional  definition  of  Z(^).   The  case  in  which  f  approaches  the 
real  axis  from  above,  which  is  of  more  physical  interest,  can  be 
studied  by  simply  taking  the  complex  conjugate  of  the  results. 

The  case  when  u  1^  1  is  of  particular  interest,  in  which  case  G 

—  1   ?    ~ 
lies  entirely  on  the  circular  arc  I,   =  R^expd'^),  for  |-^|Ctan"  {^^-   1)^ 

It  follows  from  the  symmetry  of  C,  and  the  reflection  property  Z(?)  = 

-Z*(-y*)  that  Im[l  \=0,   thus  Z    and  a are  all  real  in  this  case. 

"■  ^   '  'n  n,m      n,m 

This  property  will  be  used  in  the  following  section  in  deriving  the 

condition  for  evanescence  for  plasma  waves. 

2  -        2  - 
When  jJ^^l,  C  can  be  chosen  as  the  interval  (-= — \(~^^      •  J^^     )' 

In  closing  this  section,  one  may  point    out  an  alternative 
expression  for  the  integral  in  I^,  namely 
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III.  DISPERSION  RELATION  NEAR  CYCLOTRON  RESONANCE 

Maxwell's  equation  for  the  electric  field  can  be  written  as 

nii  (.n  y.   E)  +  6.E  =  0,  (16) 

where  n  =  kc/a3  is  the  refractive  index,  so  that  the  dispersion 
relation  for  plasma  can  be  written  as 

det(n^I  -nn  -  €)  =0.  (17) 

-  —  ~ 

Generally,  it  is  more  convenient  to  introduce  the  ordering  (9)i  together 
with  the  expansion 

E  =  E°  +'nE^  +  .  •  •  (18) 

into  Eq.(l6),  rather  than  expanding  the  determinant  (17).  One   has 
to  order  iq"^,  C    '2°  =  0,  or  ^_~^  E°  =  0,  where  it  is  assumed  that 
n  =  0(1)  0   Therefore,  E_  is  of  first  order  in  Tl  .   To  order  yT    ,  one 
has  f^'-'-^-E^  +  ef'^^.E-*-  =  0,  giving 

^i'}h\   +  ^i;,^^E°  =  0,  (19) 

and  to  the  next  order,  one  obtains 

(e(,°;-ini-n2)E°^n^,n,ei5^E°  =  0. 

(^2^-2iy.,^(_0^))^0,      -i^^O  ,  e^_-_2)E2.6(.-;)El  =  0. 


*ir-^      "II      "--     -      -II     II 


in.n,.e-tE^(e(;,)-nl)E°^^(-;)E^  =  0. 


(20) 


The  first  and  the  third  of  Eqs.(20),  together  with  (19),  form  a  system 
of  equations  for  E^,  E  ,  and  E.|  ,  the  determinant  of  the  coefficients 

yielding  the  zero  order  dispersion  relation.  [^The  second  of  Eqs.(20) 

2     1  "1 
should  be  moved  to  the  first  order  system,  since  it  contains  E  and  E  . J 
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One  has,  by  making  use  of  Eqs.(lO), 

(1  -  X  -  ni-n^XA)(2  -  X  -  nj  -  2n^ )  -  njn^  =  0,  (21) 

A  =  -*(^0,2-^0,/^0,o)-  ^22) 

The  following  comments  can  be  made  regarding  the  dispersion  relation 
(21).   First,  this  relation  differs  from  the  cold  plasma  dispersion 
relation  (evaluated  at  the  resonance)  by  the  presence  of  ain  additional 
term  in  the  first  factor,  which  contains  the  function  A.   Therefore, 
the  function  A  is  solely  responsible  for  the  generation  of  the  zero 
order  relativistic  Doppler  effects,  which  will  be  further  studied 
in  what  follows.   Second,  it  is  clear  from  Eqs.  (12)  and  (22)  that 
when  Of  ^A  1  (i.e.  n,.»'r|  )  A  approaches  zero,  which  indicates  that  the 
relativistic  broadening  becomes  negligible.   Third,  in  the  case  of 
perpendicular  propagation  (nn  =  0,  orc(->oo)  the  two  modes  represented 
by  the  two  factors  in  the  first  term  in  (21)  are  decoupled.   The 
dispersion  relation  for  the  ordinary  mode  (the  first  factor)  is 

n^  =  (1  -  X)/(l  *   XAq),  (23) 

where  A^  =  lim     A.   As  will  be  shown  in  Sec.  IV,  in  the  Maxwellian 
0    oc  ->  00 

case  this  limiting  value  coincides  with  the  result  of  Dnestrovskii  et 

2 
al.  .   The  second  factor  in  Eq.(21)   represents  the  extraordinary 

mode,  and  contains  no  broadening  effect  to  the  zero  order  theory  as 

developed  here;  thus  the  cold  plasma  dispersion  relation  applies* 

n^  =  2  -  X.  (24) 

In  order  to  study  the  broadening  of  this  branch,  the  second  order 
contributions  for  €^^   are  required.   This  will  be  done  in  Sec.  V. 
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Fourth,  as  is  demonstrated  in   See.  II,   the  imaginary  parts  of  the 
integrals  I   defined  in  (13)  are  zero  whenever  2|a^l»  indicating 
that  A  is  real  in  this  case.   Therefore,  the  region  of  evanescence 
for  waves  near  cyclotron  resonance  is  obtained  as 

|I1|/C0>1  -  in^  .  (25) 

In  the  special  case  of  perpendicular  propagation  for  Maxwellian 

plasmas,  this  result  (which  then  reduces  to  (O  <  liTl  )  has  already 

2-4 
been  pointed  out    .   It  may  be  of  some  interest  to  note  that  the 

condition  (25)  closely  resembles  a  similar  condition,  namely 

(Cl/co)  -  1  ■*■!?  >0,  in  an  earlier  work  ,  in  which  waves  propagating 

along  the  external  field  (n_L  =  0)  were  studied. 
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IV.    MAXiVELLIAN  CASE 

The  purpose  for  studying  this  important  case  is  two  fold: 
it  allows  one  to  perform  quantitative  evaluations  of  the  results 
presented  in  previous  sections,  and  it  also  enables  one  to  compare 
the  present  results  with  those  previously  published,  since,  as  was 
pointed  out  earlier,  previous  studies  invariably  assumed  this  case 
at  the  onset  of  their  formulations  by  adopting  the  continuum 
description  for  the  dielectric  tensor. 

The  relativistic  equilibrium  distribution  function  is 

F(u)  =  ikXc^)-'^[p/K^{p)]    exp(-p^),  (26) 

where  j3  =  m  c  /T   =  2r[  .   By  introducing  the  normalization  w  =  u/v  , 

2 
and  expanding  into  powers  of  ty   ,    one  obtains 

F(w)  =  K^^^    e-^^l  -  ^^2  ,  ^^2^    ,   ^^^^^^  (2^) 

Thus  1]  dependence  in  F(w)  can  be  neglected  in  the  first  order  theory. 

(2) 
One  also  observes  that  ^l^    *   which  will  be  used  in  studying  the 

broadening  of  the  perpendicular  extraordinary  mode  in  Sec.  V.,  has 

no  second  order  correction  due  to  the  relativistic  effects  associated 

with  F(w),  since  6^^/  involves  only  the  zero  order  moment,  and  <1>  =  1, 

Therefore,  one  can  use  nonrelativistic  Maxwellian  averages:  <(w  )=   i, 

<w5'^wf,"^>  =  n!(2m)!4-Vm!,  and  < w^" >  =(2n^l)  S^'Vnl .   Moreover,  Z{Z,) 

is  the  standard  plasma  dispersion  function^,  and  satisfies  Z'=-2-21:;Z, 

1 
with  Z(0)  =  iTt*  .   The  integral  in  (13)  now  can  be  performed  exactly, 

by  making  use  of  the  recurrence  relation: 

I^  =  nl^_^   -   ^C^^'Z^^)    -  ;2^^V(2nn) 

giving  n         ou  r-  t  (28) 

I^  =  -n!    I    {zf^AD  [iz(^)  -C/(2kn)] 
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Equations    (15)    now  yield 

Zq.i    =    2M|  ^|0()[z]  -¥^l^z],  (29) 

^0.2    "   20(^1^  (f  -OC^  ^io(2)    [z]    M^^^i0^2)[2,2z]   .  S0(2[^^z]  , 

where  the  notation  [fC^)]  denotes 

[FlO]  =  F»(l  -  U-2i°')*)  -   F.(^J-JL^^lM*j.         (30) 

where  the  complex  conjugate  is  taken  to  assure  that  \    approaches  the 
real  axis  from  above. 

It  is  straightforward,  but  somewhat  tedious,  to  show  that  Eqs.{29) 
reduce  to  the  nonrelativistic  limit  as  o^ ->  0,  namely 

^0.1  -  1  '?2(|)-  (3,, 

which  are  more  readily  obtained  from  Eq.(Al  )  directly,  which  states 
that  in  this  limit  Z    ^i  -»  ^wi w,,  ^ -^  Sz^   ^. 

The  case  0(  — *  od  is  also  of  considerable  importance,  since  it 

corresponds  to  the  case  of  perpendicular  propagation,  viz.,  n  = 0. 

However,  in  taking  this  limit  one  should  keep  the  quantity  ^/^  finite, 

and  consider  the  limits  ofwZ^  ^^  and  o^  Z-,  ^  ..    (cf.  Appendix  II). 

^  U , ^m        0 ,^in^l 
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Although  this  limit  can  be  carried  out  directly  using  Eqs.(29),  it  is 
more'  convenient  to  use  the  definition  (11)  directly  as  is  done  in 

Appendix  II.   The  results  can  be  expressed  in  terms  of  the  functions 

2 
F   introduced  by  Dnestrovskii  et  al.  «  as  well  as  the  functions  7" 
q  q 

3 
introduced  by  Shkarof sky-^.   One  finds  that 


^  •       (^r,  >  _  2n!  (2m)  I  „ 

^^-J^^2n.2m)--^;^7^Vm.3/2' 


(32) 


leading  to 

lim  a.   ^   -  -  — — — ■ —  r  .  j.  ^  /o  i 


T         -   2nl(2m)t  ^ 
lim  a  ^   -  -  — — * '—   F. 


lim  (o/a      )  =  -  2nl(2m-H)I(p         p       . 
^lim^((ya^^2mn^     4^  ^l     ^Vm-5/2  ^n^-7/2^- 


(33) 


The  latter  can  be  used  to  calculate  the  limiting  value  of  the  dielectric 
tensor,  cf . ,  Eq.(lO),  as  well  as  the  function  A  introduced  in  Sq.(22). 
One  obtains 

A^  =  lim  A  =  iF„/„.  (3^) 

^  oC-»-OD      ^^^ 

Equation  (3^)»  together  with  the  dispersion  relation  (23) »  coincides 

with  the  results  for  the  perpendicular  ordinary  mode  first  derived 

5  2 

by  Gershman  and  Dnestrovskii  et  al. 

For  arbitrary  values  of  c( ,  the  function  A((y,|)  can  easily  be 
calculated  by  using  Eqs.(29)  together  with  Eqs.(12)  and  (22).   Figure 
1  illustrates  the  real  and  imaginary  parts  of  A  as  a  function  of  -2^/o( 
for  various  values  of  the  parameter  ex.   The  case  U=  ^    agrees  with  the 
limiting  value  Aq  within  two  digits  of  accuracy,  leading  to  the 
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conclusion  that  the  asymptotic  expansions  in  terms  of  the  parameter 
l/b(   presented  by  Shkarofsky^,  which  are  further  studied  in  Appendix 
II,  are  useful  whenever  (x>  5«   The  cutoff  condition  (25)  is  also 
observed  in  Fig.l  as  Im(A)  vanishes  for  2^<V^1.  Both  real  and  imag- 
inary parts  of  A  monotonically  decrease  in  general,  thus  making  the 
cold  plasma  approximation  more  appropriate  for  0(  «  1. 
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V.  SOKE  APPLICATIONS 

A.  Perpendicular  Extraordinary  mode 

In  the  case  of  the  perpendicular  propagation,  the  dispersion 
relation  (1?)  can  be  factorlzed.  It  Is  more  convenient,  In  this 
case,  to  expand  directly  the  extraordinary  branch,  which  reads 

^     ^  ^  e^^  +  ie_  ^  e^_exp(-2i"//)  •  ^^^^ 

rather  than  developing  the  second  order  theory  In  the  manner  of 

Sec.  III.   Introducing  the  ordering  given  In  Eq.(9)  Into  (35).  and 

2 

keeping  terms  up  to  order  Y|  ,  one  obtains 


n2  =  2  (iO)  ,  ^2  2 


It  Is  Important  to  note  that  the  only  second  order  component  needed 

(2 ) 
for  the  extraordinary  mode  Is  6V/»  which  manifests  the  convenience 

of  the  present  coordinate  system  over  the  cartesian  system,  in  which 
every  component  appearing  in  the  dispersion  relation  must  be  evalu- 
ated to  the  second  order.   This  component  is  obtained  from  the  first 
of  Eqs. (8) ,  giving 

which  in  the  Maxwelllan  case  and  nj|  =  0  yields 

where  x^  =  -2\/6/    =  2ir|"^  (|n,|/cO  -  1 )  .  The  second  term  in  the  right-hand 

3 
side  of  (37)  is  missing  in  Shkarofsky's  "warm"  elements  .   ilowever, 

Its  contribution  to  the  dispersion  (as  will  be  sho\\Ti  later)  is  small, 

since  it  affects  only  the  Doppler  shift. 
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By  making  use  of  Eqs.(lO)  and  (33) t  and  also  approximating 
successively,  one  obtains  the  dispersion  relation  from  (36)  as 

-'  =  (2  -  X)  .  f  [%^(1  -  XF^/2)2  -  ^X(2  -  X)%/2 

+  X(|-^x2)-  ^Ai2-X)].  (38) 

The  latter  equation  has  notable  differences  compared  to  the  one 

2  k 

derived  by  Dnestrovskii  et  al.   and  used  extensively  by  Fidone  et  al.  , 

n^_p  =  (2-X)  +'^t^(2-X)V(^XFy2)-  ^39) 

The  difference  stems  from  the  neglect  of  the  second  order  contribution 

from  €.  J.,  and  also  from  the  incorrect  assumption  that  £  -  6  =  26.  e  ' 
^-  +  '  ^  XX   yy   ■*■- 

=  0.   Unlike  the  Shkarofsky  approximation,  Eq.(39)  fails  to  describe 
the  relativistic  Doppler  broadening  appropriately.   To  see  this, 
consider  the  imaginary  part  of  n  ,  which  can  be  written,  with  Eq.(A9), 
as 

Im(n2)  ='^|2(2  -  xf   ^1^(1  _  ^)2ij,(_4_)  .  |x_  ^^^^^^^j  ^       ^^^^ 

Im(n^_p)  =-rf  -^^  Im(^).  (^1) 

The  second  term  in  the  bracket  in  Eq.(40)  affects  the  damping  coef- 
ficient  only  for  larger  values  of  x, when  Im(n  )  is  already  small. 
Thus  the  major  difference  between  Eqs.(40)  and  (41)  stems  from  the 
factor  (1-2X/5)  ,  leading  to  the  conclusion  that  the  discrepancy 
becomes  increasingly  important  as  X  aproaches  its  cutoff  value  2. 
For  a  moderate  value  of  X  =  0.3  the  relative  error  is  nearly  25?^.  and 
at  the  cutoff  it  could  be  as  high  as  64^.   The  J-r  approxir::ation 
is  valid,  however,  for  tenuous  plasmas.   In  r'ig.  2,  a  comparison 
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of  the  present  result  with  the  approximation  of  Ref.2  is  presented 
for  /n  =  0.034,  ^=  (27r)l8Ghz.,  and  for  X  =  0.3  and  1.0  is  shovm.   An 
independent  calculation  provided  by  D.B.  Batchelor,   which  is  based 
on  ohkarofsky  formulation-^,  is  m  good  agreement  with  the  present 
result,  indicating  that  additional  Doppler  shift  terms  found  in 
^^•(38)  do  not  affect  the  absorption  rate. 
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B.  Nearly  Perpendicular  Ordinary  Kode 

The  appropriate  dispersion  relation  is  given  by  Eq.(21),  with 
the  definition  (22),   For  the  Kaxwellian  case,  Eqs.(29),  together 
with  (12)  are  used  to  calculate  A.  [^^  "^^^  non-Kaxwellian  case, 
Eqs.(15)  must  be  used  instead  of  (29). J  Although  the  resulting  biquad- 
ratic equation  for  n_L  can  be  solved  by  factorization,  a  simpler  approx- 
imate form  can  be  obtained  by  approximating  successively  with  respect 
to  n...   The  zero  order  solution  is 

(ni)Q  =  (1-X)/(1^AX),  (42) 

leading  to 

2 

'"i'l  =  '■^i>o[i  -  i^Z-  x)A  ]  ■  c^^) 

The  second  term  in  the  bracket  represents  the  broadening  due  to 
mixing  of  modes,  and  it  is  expected  to  generate  a  small  correction, 
since  n..  is  small. 

Numerical  evaluations  performed  for  a  Maxwellian  plasma  based 
on  Eqs.(42)  and  (^3)  are  presented  in  Fig. 3  for  various  values  of 
the  parameter  cx,  and  for  /n  =  0.03^,  (j   =(211)  l8Ghz.  ,  .(  =  0.3.   It  is 
found  that  for  oi  =  5.  both  the  real  and  imaginary  parts  of  kj_  differ 
only  negligibly  from  their  corresponding  values  for  o(=qo:   they 
appear  to  be  shifted  slightly  to  the  left.   The  coupling  effects 
are  also  negligible  in  this  case,  since  in  this  case  n..  ^i^/oc  = 
0.0068,  corresponding  to  a  propagation  angle  of  0  =  39.5  •   For 
a.=  1  (9  =  87.7  )  the  shift  is  more  pronounced,  and  the  resonance 
curves  are  flatter.   For  «  =  0.5  (6  =  85.5°)  the  broadening  is 
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further  extended  into  the  region  (ij>lfi!  i  starting  from  the  cutoff 
value  (jj  =   1.0023irif  •   The  coupling  effects  are  also  more  significant 
in  this  case,  and  the  line  width  is  more  than  doubled  compared  with 
the  case  of  perpendicular  propagation. 

The  study  of  the  nearly  perpendicular  extraordinary  mode  based 
on  the  dispersion  relation  (22)  predicts  only  broadening  due  to 
mode  coupling  effects.  Calculations  performed  indicate  that  these 
effects  are  of  the  same  (or  lower)  order  of  magnitude  as  those  obtain- 
ed in  the  second  order  analysis.   Therefore,  the  broadening  of  the 
nearly  perpendicular  extraordinary  mode  requires  the  second  order 
treatment,  which  will  not  be  given  here. 
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C.  }Iigher  Order  Harmonics 

In  order  to  study  the  relativistic  broadening  of  higher  order 
harmonics,  one  must  re-examine  the  ordering  of  the  dielectric  tensor. 
When  the  ratio  of  the  Larmor  radius  to  the  wavelength  is  of  order  ti 
(as  is  assumed  in  this  paper),  one  observes  from  Eq.(5)  that  the 
contributions  of  the  Nth  cyclotron  resonance  to  the  dielectric  tensor 
have  the  following  ordering: 

Therefore,  when  n;^2,  E  does  not  have  to  vanish.   For  N]^  3,  the  zero 
order  components  do  not  contain  the  resonance  effects,  thus  the  cold 
plasma  dispersion  relation  evaluated  atW=Nfn.|  prevails,  viz., 

-  nnnid  -^   -  n^)   =0.  {^k) 

'      N^-1 

For  N=2,  however,  the  dielectric  tensor  contains  a  resonance  contri- 
bution through  the  component 

e__  =  1  -  ix  *^iXni  a-^^Q  +  0(r|),  (^5) 

which  modifies  the  cold  plasma  dispersion  relation,  yielding 

Dg  +  iXni  a^^Q  V.^,    =  0,  (46) 

where  M,,  is  the  cof actor  determinant  of  the  (--)  element,  viz.  , 

The  function  a.  „,  which  is  defined  in  Eq.(ll),  can  be  written,  by 
1 ,  u 

using  Sqs.(Al)  and  (A2)  as 
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where 

The  functions  Z    can  be  calculated  from  Eqs.(15)i  or,  in  the 
n,m  -1   X  ^/  f    t 

r.laxwellian  case,  from  Eqs.(29). 

2 
The  dispersion  relation  (46)  is  a  cubic  polynomial  in  nj_,  due 

to  the  presence  of  the  relativistic  term  a.    ^.      However,  in  the  case 

of  perpendicular  propagation  it  can  be  factorized.   The  first  factor 

leads  to  the  cold  plasma  dispersion  relation  for  the  ordinary  mode 

evaluated  at  the  resonance,  i.e.,  • 

n^  =  1  -  XA,  (50) 

indicating  that,  unlike  the  case  of  N  =  1,  the  second  harmonic  ordinary 
mode  contains  no  zero  order  broadening.   The  extraordinary  mode, 
which  is  associated  with  the  second  factor,  is  obtained  from  the 
following  dispersion  relation: 

(l-X/6)(l-x/2)  -n2(l-X/3)  ^n^X  (l-xA-in^)  H^^^-   0,  (51) 

where  a.  ^  is  the  large  o(    limit  of  a-  q,  which  in  the  ^':axwellian  case 
gives  a?  Q  =  -ZF^/p,  cf.,  Eq.(Al4).   Therefore,  for  N  =  2,  the  absorp- 
tion of  the  extraordinary  mode  is  more  significant  than  that  of  the 
ordinary  mode  when  n..  =  0  (or  nearly  zero).   Apart  from  notation, 
Eq.(51)  in  the  I-'iaxwellian  case  is  equivalent  to  those  given  in  Refs. 
2  and  3,  and  to  that  recently  studied  by  Bornatici  et  al.    in  connec- 
tion with  certain  mode  conversion  schemes. 
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The  two  factors  in  the  first  term  of  the  dispersion  relation  (51) 

correspond  to  the  cold  plasma  right  and  left  hand  cutoff  frequencies 

2 
at  the  resonance,  and  the  coefficient  of  n   in  the  second  term  leads 

to  the  upper  hybrid  resonance  frequency.  The  last  term  generates  the 

relativistic  broadening  effects. 

When  n  .  is  not  necessarily  zero,  Sq. (^6)  can  be  written  as 

A^y^  +  A2y^  +  A^y  +  A^  =  0,  (52) 

where  y  =  nj_,  and  with  the  shorthand  X  =  1  -  X/m,  the  coefficients  are 


m 


A^  =  i  X  a^^Q, 

A2  =  X^  -  A^(X^  +  2X^  -  nf,). 

A^  =  -X^X2  -  X^(X^  -  n^)  ^  X^4  -  A^  2X^  (X^  -  n^), 

In  Fig.  4  the  function  a..    ^   is  illustrated  for  the  Maxwellian  case 
for  the  values  of  «  =  5i  1»  and  0.5.   The  imaginary  part  of  a.  q  is 
zero  when  2fp0(^  1,  thus  the  region  of  evanescence  for  the  second 
harmonic  is 

2{al/(v^    1  -  i-n^i-  (53) 

For  c(=5i  a.  ^  is  very  close  to  its  large  a  limit  -2F„/2,  and  for 
oc  «  1,  both  real  and  imaginary  parts  of  a.  ^  are  much  less  than  1 
in  magnitude,  hence  the  broadening  is  negligible  in  the  latter  case. 

Figure  5  illustrates  the  variation  of  the  real  and  imaginary  parts 
of  kj_  for  the  extraordinary  mode  as  a  function  of  2|ri.l/6;,  for  the 
r.:axwellian  case  with  ;:  =  O.3,  v  /c  =  0.03^,  and  W  =  27T  18  Ghz. 
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The  general  behavior  of  these  curves  is  similar  to  those  given  in 
Fig.  3  ^or  the  fundamental  ordinary  mode.   The  c<  =  5  case  coincides  with 
the  perpendicular  propagation  case  within  at  least  three  decimal 
digits.   For  smaller  values  of  o< ,  the  maximum  value  of  the  imaginary- 
part  of  kj_  decreases,  and  the  line  width  increases. 

The  second  harmonic  ordinary  branch  has  a  negligible  imaginary 
part  and  does  not  vary  appreciably  from  the  cold  plasma  value  given 
in  Eq.  (50)  •   It  is  important  to  note  that  the  Re  kj_  curves  for  the 
ordinary  and  the  extraordinary  modes  intersect  at  the  two  values  of  the 
wave  frequency.   The  one  corresponding  to  the  smaller  frequency,  (or  the 
larger  value  of  2 IH. /cj) »  appears  to  be  promising  for  a  possible  mode 
conversion  between  these  branches. 

The  third  root  of  Eq.  (52)  leads  to  a  branch  which  is  nearly 
electrostatic  and  heavily  absorbed  in  the  plasma. 
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VI.  SUMMARY  AND  DISCUSSIONS 


In  the  present  paper,  the  relativistic  broadening  of  nearly 
perpendicular  waves  near  cyclotron  resonance  is  investigated  for 
a  warm  (v  /c^l),  isotropic,  but  not  necessarily  Maxwellian, 
plasma.   The  underlying  assumptions  are  that  a)=0(CiL)  =0(12)  =0(kc), 
so  that  £  1=  0(c  /v  )  and  n..  =  0(v  /c).   When  these  conditions 
are  satisfied,  one  also  has  j^-.  |  =  0(c/Vg),  and  the  zero  order 
dispersion  relation  is  no  longer  the  cold  plasma  dispersion 
relation  near  the  resonance  as  was  found  in  the  nonrelativistic 
case.'^  The  formulation  is  uniform  in  the  parameter  v^/(cn^V  so  that 
all  values  of  n,.  [provided  that  n  =  0(1)]  can  be  studied.   However, 
the  relativistic  effects  are  unimportant  when  this  parameter  is 
much  less  than  one.   The  formalism  used  is  also  suitable  in  geo- 
metrical optics  methods  to  study  the  effects  of  weak  inhomogeneities, 

even  though  the  magnitude  of  these  effects  may  be  altered  due  to 

7 
the  ordering  assumed. 

The  basic  advantage  of  the  formulation  given  here  is  to  express 
the  dielectric  tensor  in  terms  of  a  set  of  functions  defined 
similar  to  the  plasma  dispersion  function  of  the  nonrelativistic 
theory.   Various  recurrence  relations  and  reduction  formulae 
obtained  help  one  to  evaluate  these  functions  with  relative  ease; 
particularly  so  for  the  Maxwellian  case,  in  which  case  the  results 
are  expressed  without  any  integrations  or  infinite  series  involving 
the  plasma  dispersion  function.   The  latter  property  enables  one 
to  calculate  efficiently  the  broadening  of  the  plasma  spectrum  in 
some  range  of  physical  space. 
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The  formalism  is  applied  to  the  fundamental  extraordinary  mode, 

(which  requires  the  second  order  theory,)  for  n.j  =0,  and  the  correction 

2 
terms  for  Im(n  )  are  given,  which  generalize  the  previously  published 

results  for  the  tenuous  plasmas.   The  study  of  the  fundamental 

ordinary  mode  shows  that  the  relativistic  broadening  is  very  sensitive 

to  the  angle  of  propagation,  and  it  is  restricted  to  the  region  within 

a  few  degrees  from  the  perpendicular  propagation. 

2N-4 
The  broadening  of  the  higher  order  harmonics  is  of  order  y) 

for  the  Nth  harmonic.   The  zero  order  theory  for  the  second  harmonic 

shows  that  the  ordinary  mode  is  nearly  undamped  and  the  extraordinary 

mode  has  a  damping  rate  larger  than  the  one  for  the  fundamental 

ordinary  mode.   The  possibility  of  the  mode  conversion  between  these 

branches  are  also  pointed  out.   An  analysis  for  the  latter  problem 

requires  the  inclusion  of  the  spatial  inhomogeneities,  which  will 

be  studied  elsewhere. 

The  region  of  evanescence  for  the  Nth  harmonic  is 

N  lnl/60  <  1  -  inf,.  (5^) 

It  is  of  interest  to  note  that  this  condition  is  independent  of  the 
plasma  frequency,  thus  it  may  be  used  as  a  diagnostic  tool  in  weakly 
inhomogeneous  plasmas. 
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APPENDIX  I:  PROPERTIES  OF  Z     FUNCTIONS 

n,m 

The  functions  Z    are  defined  in  Eq.(ll)  with  D  given  in  (3) 

are  needed  only  for  even  values  of  n.   The  following  recurrence 

relation  is  obtained  directly  from  the  definition 

V2.m=^2/o()(Z^^^,,-^Z^^^-  <wM»-Z^^^,2.  (Al) 

which  enables  one  to  express  these  functions  in  terms  of  Z„ 

^  0,m 

Furthermore,  when  F  is  isotropic,  as  is  assumed  in  this  paper,  the 
derivatives  of  these  functions  can  be  calculated  using 

dZ   /d|  =  Z'   =  mZ     ,  -nZ   „   ^,,  (n>0),  (A2) 

n,m   *•    n,m    n,m-l    n-2,m-^l'  ^  "  '  '  ^   ' 

which  is  obtained  by  integrating  by  parts  after  replacing  F  in  (11) 
by  its  derivative. 

Since  D  =  -|«(w  -  ^^)  (w  -  X,_)  ,    where 

X_  =  (i/o^)[e  ^  (e^  -2fo()^]  .  (A3) 


u  being  the  cosine  of  the  polar  angle  in  the  velocity  space,  one 
obtains,  with  the  use  of  the  identity 


k-1  k-^1         ^ 

w 


w  -  a       w  -  a        ■  _^ 

J    ^  r-  y2n+m+2  y2n+m+2 

2n,m             Qf-'-lq            ^  CC                  [w-'C          w-^ 

2n^-n  -1 

where  q  =  ( l/oc)  (^^  -  2f  o()  ^  =  i(C^  -  ^_)  .   Since  C  ^  "B)  =  "^^(f )  '  after 

extending  F  as  F(-w)  =F(w)  and  replacing  y ->  -u  and  w-»-w,  one 
obtains 
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^  2n+m  j 

where  the  function  Z(C^)  is  as  defined  in  Eq.(l4).   Equation  (A5)  can 
be  further  reduced  by  observing  that  the  lower  half  of  the  y  integra- 
tion in  the  first  term  of  (A5)  is  (-1)"^  "^  times  that  of  the  upper 
half,  hence  j+m  must  be  even.   Moreover,  since  dyi/q  = -<Xd^ +/^+,  with 
0^^+  -  2^iC^+  *-  2\   =0,  the  two  terms  in  the  second  integral  can  be 
combined  into  one,  by  introducing  the  variable  Z,   as 

C^'t;  =  ^  +  {/  -2'^o()^,  (A6) 

and  the  contour  C  defined  as  the  image  of  -l^u</l  under  the  mapping 
(A6).   One  obtains,  with  r=rm/2]  being  the  closest  integer  less  than 
or  equal  to  m/2, 

p  0  J  -0 

-  2  Jd^[c2-(|0((2,^)2jn(.p^^2,|)m^2(^).       (A?) 

Equations  (15)  are  obtained  directly  from  (A?).   Since  the  factors 
which  multiply  l^Zil^)    in  the  integrand  in  (A?)  are  polynomials  in  K   , 
the  dependence  of  the  second  term  on  the  functions  I   introduced 
in  (13)  is  evident. 
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APPENDIX  II.  CONNECTION  WITH  F   and  ^     FUNCTIONS 

q    q 

In  this  Appendix  the  relations  between  Z    functions  and  the 

ii  I  In 

F  functions  introduced  by  Dnestrovskii  et  al.  is  given.   It  is  also 
shown  how  they  are  connected  to  Shkarof sky's  7  functions,  so  that 
the  results  presented  in  this  paper  can  be  compared  efficiently  with 
those  previously  known. 

Consider  the  functions,  with  F(w)  =  7l"-^  exp(-w  ) , 

^n+3/2    2n+l  !!  V        i    2      2xm  *  ^^^' 

J'  I         i  (w  -X  ) 

These  functions  have  the  following  properties «   For  m  =  1, 
F^j>]  =  2  *2xZ(x)  =  -Z'  (x)  ,  with  x  =  -(x^)^, 


F^J)  =  (l^xM^))/.. 


(A9) 


Since  these  are  the  recursive  properties  which  generate  the  Dnestrovskii 
functions,  one  concludes  that 


q 

For  m>  1, 


F^^^  =  F  .  (AID) 


F^';)  =  [(m-l)F^^-l^+  x2F^"^)]/q. 


__d  p(m)  ,  p(mn)  (^11) 

d(x2)  q      q 

p(m+l)  _  p(in)  _  p(in) 
q+1     q      q+1 

In  order  to  relate  these  functions  to  Z^  ^^    consider  the  asymptotic 
expansion  of  the  latter  in  l/«  which  is  obtained  by  expanding  1/D 
into  power  series,  giving 


z 

n,m 
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^  2k-m+l  f  ^   w^w2^F(w) 

^^^  ,  2    2,2k-mn  •  ^^12) 

(w  -  X  ) 


uu 


k=r 


where  r=[(m+l)/2]  ,  and  x^  = -2|/o(  =  2y|"^  (|Xll /(J  -  1)  • 

One  obtains,  after  some  straightforward  manipulations,  and  angular 
integrations. 


—       (2k)I 
^2n. 


'^  2"^-^  k=r  CX^^-^^^^k!  C2k-m)!   ^^^^3/2  '  '^^^^ 


The  large  oc  limits  given  in  Eqs.(32)  follow  directly  from  (AI3), 
with  the  third  of  Eqs.(All).   Substitution  of  (AI3)  into  the  defini- 
tion of  a    in  (11)  yields  the  result 
n,m        "^ 

-   nJ_  y  ^^^^'  p(2k-m^l)  .,.K^ 

^n,m--  ^m-l  ^^   ^2k-m   ,,  ^^k-m)  !   ^^^^^5/2   ' 

from  which  the  limits  in  Eqs.(33)  follow.   When  this  result  is 
substituted  into  Eqs.(lO),  one  finds  that  Shkarofsky's  "warm"  elements 
are  obtained  to  the  order  one  in  in   and  order  -2  in  p<. 

In  order  to  establish  the  relationship  with  Shkarof sky ' s  T 
function  (or  W  function  used  in  Ref.^)  consider  the  simplification  (15) 
given  in  Ref.3 

^^  2  io   (i_i,)q^"H^^  0(2(l-it)^ 

=  exp(-l/a2)  ^  -^ —  F;5ihx2^  l/a2).  (AI5) 

j=o  cx  -^  jl   ^  J 

In  the  latter  equation,  the  function  F^  ^    is  regarded  as  a  function 

2 

of  X  .   Expanding  this  function,  as  well  as  the  multiplicative  expo- 
nential function,  into  powers  of  l/a,  and  rearranging  the  terms 
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one  obtains  the  asymptotic  series 

CD 

:k+m) 


K""^ '  L  ^^  ^iVi  ■■■' .  (A16) 


^y    =  ?  .  and   T:'-'!'  ^-^F'"'.  (A17) 


where 

^q      ■'q  •  -^'-   ^q       ^^^2^  ^q 

In  deriving  (Al6),  the  following  relation  is  used,  which  can  be 


established  by  induction, 

k 

(A18) 


k  p 

p(2k*m)  ^  y  /.ik-p  fk\y  (p\  p(p+m-j) 


p=0        ^i^^j^O^J^  ^^J 

The  functions  a    for  m  =  0,l,  and  2,  which  appear  in  Eq.(lO),  can 
be  expressed  in  terms  of  Shkarofsky  functions  as 

^n,l  =  -^2nl/a)7^^2)/2.  (A19) 

Equations  (A19)  can  be  further  reduced  by  making  use  of  the  recurrence 

relation 

^(m-^1)  ,  ^m)  _  c^m)  ^  (^20^ 

-'  q  +1    -'  q     ^  q+1 

The  function  A  may  also  be  written  as 

An  evaluation  of  A  based  on  the  latter,  together  with  (A15) .  is 
certainly  much  less  efficient  than  the  method  presented  in  this 
work.   Keeping  terms  up  to  order  l/b^  in  (A21),  fusing  (Al6)  ,J  and  substi- 
tuting into  (^2),  one  gets  the  dispersion  relation  given  in  Ref.  3« 
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FIGURE  CAPTIONS 

Fig.  1     Real  and  imaginary  parts  of  A(f,0()for  the  Kaxwellian 
case  as  a  function  of  x  -  -2|/X  =  2^^  (  lo/col  -  1)  ,  for 
CC=  5.0,  1.0,  and  0.5. 

Fig.  2     A  comparison  of  Dnestrovskii  approximation  with  the 
present  analysis:  Im(k_L)  (cm"  )  vs.\SX/c^\     for  the 
perpendicular  extraordinary  mode,  and  for  X=0.3,  1.0, 
1^=0.03^,  0)=  (27i)l8Ghz.  . 

Fig.  3     Real  and  imaginary  parts  of  kj_  (cm"  )  as  a  function  of 
(iV60|  for  the  nearly  perpendicular  ordinary  mode,  for 
X  =  0.3.  '^=  0.03^1-,  60=  (2T()  18GHz.  ,  and  for  the  values 
of  o<^  =  0.5,  1.0,  5.0.  (Kaxwellian  case) 

Fig.  4     Real  and  imaginary  parts  of  a.  ^  for  the  Maxwellian 

case  as  a  function  of  x^  =  -2^2/«- ""  2r|^(2ll^/u)  -  1),  for 
C'.=  5.0,  1.0,  and  0.5. 

Fig.  5     Real  and  imaginary  parts  of  k_L  (cm"  )  as  a  function  of 
2  lill/cJ  for  the  second  harmonic  extraordinary  mode,  for 
X  =  0.3,  00=0.034,  CJ  =    (27l)l8Ghz.,  and  for  the  values 
of  oc  =  0.5,  1.0,  0.5.  (riaxwellian  case) 
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